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The features of alternative approach of non-equilibrium evolution thermodynamics are shown on 
the example of theory of vacancies by opposed to the classic prototype of Landau theory. On this 
foundation a strict theory of the melting of metals, based on development of Frenkel ideas about 
the vacancy mechanism of such phenomena, is considered. The phenomenon of melting is able to be 
described as a discontinuous phase transition, while the traditional Frenkel's solution in the region 
of low-concentration of vacancies can describe such transition only as continuous one. The problem 
of limiting transition of shear modulus to zero values in the liquid state, as well as the problem of 
the influence of extended state of vacancies on their mobility, is discussed. 
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I. INTRODUCTION 

The solids have a large number of variants of crys- 
tallographic structures and violations of these structures 
with such defects as vacancies, dislocations, grain bound- 
aries etc. Under different severe external influences, such 
as radiation, severe plastic deformation and others they 
change the internal crystallographic and defect struc- 
tures. To describe the changes the techniques of the the- 
ory of phase transitions is used, where the concept of 
order parameter or internal state variable is introduced. 

The rudiments of the concept of internal state vari- 
able were offered yet Duhem in 1903 [^ (see also more 
detailed review in Ref. [2J). In 1928 the idea was rean- 
imated by Herzfeld and Rice for the account of internal 
structure evolution of polyatomic gas when they stud- 
ied the dispersion and attenuation of sound Q. It got 
development in other numerous researches of that time 
[^-i2]. Supplemented with elements of rational mechan- 
ics, the concept took a completed and closed form in the 
works by Coleman and Gurtin 'S, ^ . In the modern sci- 
ence, the line is developed too lO-ldJ. In the works by 
Landau the concept was developed anew, through intro- 
ducing the internal state variable in the form of the or- 
der parameter, the theory of phase transitions was built 
by him at a high-quality level [1^, [la| . The Landau ap- 
proach has resulted in further development in the theories 
of the phase fields p7H25| . In addition, the mesoscopic 
non-equilibrium thermodynamics approach was develop- 
ing to fit the description of soft matter |;26n2S | . Studies 
of the problem gave birth to a plenty of variants of kinetic 
equations, describing the evolution of internal structure 
during non-equilibrium processes within the framework 
of different approaches (see [S^ and review there). 

Earlier, an alternative approach of non-equilibrium 
evolutional thermodynamics (NEET) proposed to de- 
scribe the processes in solids with structural defects (301 - 
[33| . In Ref. [30| the problem was subdivided into thermal 
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and structural parts. In the first part, the molecular dy- 
namics simulation was used to represent a thermo-motion 
as a superposition of equilibrium and non-equilibrium 
contributions. The non-equilibrium constituent was the 
acoustic emission, arising from dynamic transition phe- 
nomena at a moment of defect (dislocations) formation. 
The acoustic emission has not been considered such con- 
text before. There it was considered not only as an at- 
tenuating wave process as a consequence of both scat- 
ter on internal heterogeneities and nonlinear pumping to 
other frequency regions of the spectrum, but as a part 
of thermodynamic process of the internal energy trans- 
formation. The resulting were kinetic (evolution) equa- 
tions for the production of non-equilibrium entropy and 
its sinking into an equilibrium thermal subsystem. 

The structural part was demonstrated on the exam- 
ple of a solid with vacancies, and the evolution equa- 
tions were deduced in terms of the internal energy. In 
Ref. [3^] this part was justified with the attraction of 
the elements of statistical physics, namely, the probabil- 
ity distribution function (PDF). There are two levels of 
rigor in the description of non-equilibrium processes. At 
the first level the PDF, which does not change in time 
is taken. Maxima of this function determine the equilib- 
rium (most probability) states and the system tends to 
one of them starting from an arbitrary non-equilibrium 
one. The second level takes into account the possibil- 
ity of probability evolution in time (kinetic equations of 
Boltzmann, Fokker-Plank etc) . But already at first level 
of rigor one can obtain new additional relationships. 

There are few approaches to description of the melting 
of solids. The most simple is based on the Lindemann 
criterion, which determines a melting-point by the am- 
plitude of thermal vibrations. That is, when the ampli- 
tude becomes of the order of 10 percent of interatomic 
distance, the crystal lattice loses stability and collapses 
(melts) [3^,[3a]- This approach is not, however, universal 
type and for different crystals of BCC-, FCC- and HCP- 
type, the Lindemann constant is different. In addition, 
the Lindemann theory is one-phase, it does not deter- 
mine the free energy for the liquid state [33| . The second 



approach is based on the mechanical Born criterion of 
stability, a solid melts when the shear modulus becomes 
zero [33, |3^- However, the Born model is also one- phase 
theory which does not contain direct description of the 
liquid phase. The third approach, issuing from Frenkcl 
and Eyring works, is based on order-disorder transitions 
with the participation of vacancies [40l - |42| . Unlike the 
former approaches this theory is biphase, however, up to 
now the discontinuous nature of transition between the 
solid and liquid state was not taken into account. 

In this paper, within the framework of the developed 
NEET approach the vacancy theory of the melting of 
solids taking into account its discontinuous nature is con- 
sidered. In part II, formulation of NEET is given for a 
solid with vacancies and interstitial atoms. For a homo- 
geneous problem the PDF is obtained and a new form 
of kinetic equations symmetric in view of the use of free 
and internal energy is deduced. In the same part the va- 
cancy theory of the melting of metals is considered from 
positions of Frenkel's concept. In part III, the discus- 
sion of the obtained results is given. The problem of the 
shear modulus tending to zero during transition to liq- 
uid media is discussed. Part V contains the summarizing 
conclusions. 



II. SOLID WITH VACANCIES AND 
INTERSTITIAL ATOMS 



configurational entropy Sd is also a one- valued function 
of the number of interstitial atoms, which does not also 
depend on the energy of interstitial atoms and temper- 
ature. Provided that the microstates of vacancies and 
interstitial atoms are statistically independent, the total 
thermodynamic probability W of the microstates is 

W = WvW, (5) 

and the total configurational entropy Sc equals to the 
sum of entropies of subsystems 

Sc = Scv + Sci- (6) 

Total probability of the state containing ny vacancies 
and rii interstitial atoms will contain a limiting exponen- 
tial multi plie r, probability of Gibbs for this microscopic 
state (13,111 
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(7) 



where Z is the normalizing statistical sum, U{nv,ni) is 
the internal energy taking into account the existence of 
vacancies and interstitial atoms, T is the temperature. In 
conditions of the total degeneration the internal energy 
is a linear function of the defect number |40| 



C/ = C/q + uvonv + UiqUi 



(8) 



A. Probability distribution function 

The expressions for the thermodynamic probability 
Wv and the configurational entropy Scv assuming the 
maximal degeneration of the microstates were obtained 
by Boltzmann for a solid with vacancies |4(j | 



Wv = 



{N + nv)\ 



N\nv\ ' 
ScV = ks'inWv, 



(1) 

(2) 



where A'^ is a total number of atoms in the crystal, ny 
is a number of vacancies, ks is the Boltzmann constant. 
Note that the configurational entropy Scv is a one- valued 
function of the number of vacancies and it does not de- 
pend on their energy (and on temperature as well). 

The expression for the thermodynamic probability Wi 
and the configurational entropy Sd for interstitial atoms 
can be quite similarly written down for the maximal de- 
generation in microstates 



W,= 



iN, + n,y. 



Sc, = kBlnW,, 



(3) 
(4) 



where Ni is a total number of equilibrium positions for 
interstitial atoms in a crystal (local minima of the po- 
tential energy), rij is number of interstitial atoms. For 
different symmetries of crystals the number of such posi- 
tions can be different and multiple to N. Note that the 



where Uq is the internal energy with no contribution of 
defects taken into account, uyo, Uio are average ener- 
gies of a vacancy and an interstitial atom, which in this 
case are constant values for all possible configurations of 
atoms of the system. 

The condition of total degeneration is practically ex- 
actly satisfied for the low numbers (concentrations) of 
vacancies and interstitial atoms, when the interaction 
between defects can be neglected. At the same time, 
this condition is not fulfilled for those configurations for 
which vacancies are close to each other or merge at all 
(bi- vacancies, triple vacancies, vacancy pores) and at a 
high vacancy concentration. The same refers to the in- 
terstitial atoms as well. In these cases it is necessary 
to take into account the removal of condition of total 
degeneracy due to the interaction of point defects. As 
known, the coupling interaction reduces the total energy 
of a system that reduces the effect of action of the lim- 
iting exponential multiplier in Eq. [3] and results in a 
possible appearance of long-living configurations, which, 
thus, can compete with more numerous configurations of 
higher energy. 

The probability distribution function with total (inter- 
nal) energy Ei is [4^ 



/(i?0 = ^exp(---). 



(9) 



where w{Ei) is the distribution of states in energy Ei or 
the number of microstates (configurations) with energy 



El, Z IS the statsum over all energy states of the system 
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The states are numbered in the order of energy growth 
Ei+i > El. Distribution w{Ei) depends on the number 
of vacancies and interstitial atoms, and on their ratio, as 
well as on symmetry of their location (ordering) . On this 
account, it is impossible to fix one-to-one correspondence 
between internal energy and number or concentration of 
defects. Finding of function w{Ei) is a hard combination 
problem, determination of power spectrum of Ei as a 
set of acceptable energies for I state, is also a stubborn 
problem. However, such one-to-one correspondence can 
be obtained for the average value of internal energy U 
over all of the states for the fixed number of vacancies 
ny and interstitial atoms Uj. Using expression for PDF 
©J it is possible to write down 



1=1 ^ 



(11) 



Here, unlike the case of Eq. ([8|) , the internal energy is 
not the linear function of ny and n^, but is the function 
of general type. The increase of fraction of the sym- 
metric low-energy states diminishes the average value of 
the internal energy. The transition to the states, where 
high-concentration vacancies and interstitial atoms can 
pass to higher-energy extended (and mobiler) states [4Q- 
|48| . results in internal energy growth. To reflect all the 
properties of the internal energy it is presented in the 
form of a polynomial with alternating signs 



U = Ua + 
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The first sum describes a contribution from every sub- 
system of point defects separately; the second one de- 
scribes a contribution from the interaction between va- 
cancies and interstitial atoms. Sign "minus" at the sec- 
ond sum is because the interaction between vacancies and 
interstitial atoms is, on the whole, of the attracting na- 
ture. Ignoring the degeneracy within one set of the states 
with the same number of vacancies ny and the number 
of interstitial atoms n^ , for such average internal energy 
we can write down the effective probability distribution 
function in the form of Eq. [7] with determinations ([T]) , 
([3]) and ([5]), in which the internal energy is set not by 
Eq. ([5]) , but by a more general expression p^ . Bringing 
the variables independent of the number of defects ny 
and rii into the "inessential" constant 1/ Z we can write 
expression ([7]) in the form 



N Ni 

f{ny,n^) = —'[[{nv + i)Y[{nt+j)exp{ 
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k^ 
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Note that the Boltzmann thermodynamic probability 
for vacancies and interstitial atoms is built within the 
concept of statistical independence of forming these sub- 
systems (positions of interstitial atoms do not coincide 
with positions of lattice sites). At the same time, in the 
total distribution function these systems are not indepen- 
dent because of the mixed terms in the internal energy 
(see the last sum in p^ ). 

Values of ny and n^, at which PDF has extreme val- 
ues (equilibrium states) , can be found from the following 
transcendent equations [34] 
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where 
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are energies or chemical potentials of a vacancy and an in- 
terstitial atom. Eqs. ([TS]). p7)) are the equations of state 
for a general non-equilibrium case. Eqs (fT4|) . (flS)) mesh 
due to the mixed terms in the expansion of the internal 
energy in (J12p . determining the interaction between the 
subsystems of vacancies and interstitial atoms. 

The first two terms in Eqs dTi]) . (|T5|) are sums S of 
slowly divergent harmonic series. This part of the equa- 
tions depends only on system size and doesn't depend on 
material parameters. It is the fundamental part of Eqs. 
([T4| , (fTSJ) , decreasing with the growing parameter ny or 
rii, accordingly (curves 1 of Fig. 1, 2, 3). To calculate 
the fundamental curve S we took N — Ni = 2000. The 
last terms in Eqs. (ITi|) . (|15l) depend on the parameters 
of material through the coefficients of Wfc, and they are 
materiel parts (M) of Eqs. p^ . ([15]). As seen, the posi- 
tions of Eqs. p4p . P^ roots coincide with the maxima 
of PDF, calculated directly by formulas dH]), P^ . 

In the case of low concentration of vacancies and in- 
terstitial atoms Eqs (dU), p^ can be solved analytically 
and it gives the well-known Frenkel solution p(]| . 



B. Equilibrium states 

It is known that the interstitial atoms have higher mo- 
bility as compared to vacancies, therefore they come to 
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FIG. 1. Solution of Eq. (|14ll for a constant energy of vacancy. 
The fundamental curve S{1), the material curve (2) (propor- 
tional to energy of vacancy uv) and PDF (3): uvo = ^(s-V). 
Other parameters equal zero. T = 2400; 3000; 4000(7^). The 
right-hand scale is for PDF. Hereafter the graphs are con- 
structed as a function of vacancies concentration in percent 
n = lOOn/N. 



FIG. 2. Solution of Eq. p4|) for linear dependence of va- 
cancy energy versus the vacancy concentration. The fun- 
damental curve S{1), the material curve (2) and PDF (3): 
Uvo = 1(6^), uvi = 0.0018(61^). Other parameters equal 
zero. T = 2400; 3000; 4000(ii:). The right-hand scale for 
PDF. 



the equilibrium state first (adiabatic limit). For the times 
larger than the adiabatic limit the contribution of inter- 
stitial atoms can be ignored, as a result we get a pure 
vacancy problem [Sj. Consider this problem in varying 
degree of approximation. In the linear with respect to 
the number of vacancies approach for the internal energy 
or for the constant vacancy energy there is only one solu- 
tion of Eq. (fT4|) that is the well-known Frenkel solution 

M (Fig- m 



nve 



Nexp{- 



knT 



(18) 



where nve is the equilibrium value of the number of va- 
cancies. 

At lowering the energy of vacancy or at the increase 
of temperature the solution is continuously displaced to 
region of higher number (or concentration) of vacancies. 
We can formally get any value of vacancy concentration, 
but, it is clear, that too high vacancy concentrations of 
about 100 percent lose physical sense and the theory be- 
comes useless. However Frenkel proposed to consider not 
too high vacancy concentrations of about 10 persent as 
transition to the liquid state. At the same time, the 
solution proposed by him (in approximation of the inde- 
pendence of the vacancy energy of the vacancy concentra- 
tion) can not explain or describe the process of melting of 
a solid as a jump-like first-order phase transition. Tran- 
sition to the region of high concentration of vacancies is 
continuous during temperature growth. 

In the quadratic with respect to the number of vacan- 
cies approximation for the internal energy or in the linear 
approximation for vacancy energy Eq. (|14p has already 
two solutions in the region of high energy of vacancy or in 
the region of low temperatures (Fig. [5]). One of them in 



the region of a lower concentration of vacancies, actually, 
coincides with the Frenkel solution, but it is somewhat 
shifted due to the interaction of vacancies. The second 
solution in the region of a higher concentration of vacan- 
cies describes the equilibrium or stationary solution of 
the problem as well; however, it corresponds to a mini- 
mum of PDF and is unsteady. Probability of this state 
is lower not only as compared to the stable stationary 
state but also to any non-equilibrium state. In addition, 
for the low energy of vacancy or for high temperatures 
we enter in a region in which the solutions of Eq. (jl4[) 
are absent at all. The last circumstance testifies unsuit- 
ability of the approximation in this region, so a higher 
approximation is needed. 

In the cubic with respect to the number of vacancies 
approximation for the internal energy or in the quadratic 
approximation for the energy of vacancy, Eq. (1141) can 
already has one or three solutions (Fig. |3]). In region of 
high values of the energy of vacancy or low temperatures 
Eq. p4)) has one solution, which is a Frenkel solution 
modified due to the nonlinear contributions. At lower- 
ing the energy of vacancy or with temperature growth 
the equation has three solutions, one of which (the left- 
hand) is the modified Frenkel solution, the second one 
(intermediate) is unsteady, and the third one (the right- 
hand) in the region of high concentration of vacancies can 
be treated as one proper to the liquid state of the matter. 
In this case, we have an equilibrium coexistence of solid 
and liquid phases of the matter. For a still greater in- 
crease of temperature we pass to the region with only one 
solution, which corresponds to clear-melted matter. The 
transition from the solid state to liquid one is a jump-like 
first-order phase transition. In order to talk about the 
transition just to the liquid state (not for example to the 
amorphous state), we need not only the high concentra- 




FIG. 3. Solution of Eq. (|14p for the quadratic dependence of 
vacancy energy on the vacancy concentration. The fundamen- 
tal curve 5(1), the material curve (2) and PDF (3): uvo ~ 
0.125(eF), uvi = 0.6 ■ lO-^(eV), uv2 = 0.125 • lO-'^(eV). 
T = 260;286;300;330(ii'). The case of T = 286^" is shown 
by dotted line. The right-hand scale is for PDF. 



tion of vacancies, but also the shear modulus of material 
going to zero. This issue will be discussed in the next 
part in detail. 

Influence of the interstitial atoms can be qualitatively 
estimated understanding that they have higher energy as 
compared to vacancies, and therefore in the equilibrium 
state have more low concentration. As, with (1131) and 
(J17p . the sign "minus" is at the constant un, then, due 
to cross effects, this results in diminishing of the energy of 
vacancy uv, and to some displacement of the roots of the 
equation (|14p to the region of a higher number (concen- 
tration) of vacancies. For the case of interstitial atoms to 
be directly participating in the processes of melting the 
equation (fTS]) should have one more solution. As their 
energy is higher, their material curve will be consider- 
ably higher than that for vacancies. Consequently, while 
the equation P^ has two solutions, the equation ((T5]) 
has only one solution, and the interstitial atoms will not 
directly and independently participate in the processes 
of melting, but will only modify the participation of va- 
cancies a little. 



C. Set of non-equilibrium thermodynamic 
potentials 

The free energy F ~ U — TS, where S is the entropy, 
was introduced for solution of quasi-equilibrium thermal 
problems in the defect-free condensed matter, as one of 
types of thermodynamics potential, along with the in- 
ternal energy, enthalpy and Gibbs energy. In certain 
external conditions it, as well as other thermodynamic 
potentials for other conditions, possessed extreme prop- 
erties as opposed to internal non-equilibrium processes. 
For imperfect solids, in particular solids with vacancies. 



according to the hypothesis of Boltzmann, the free (con- 
figurational) energy was introduced by analogy 

F,v = U^ TScv, (19) 

where configurational entropy Scv is set by Eqs ([IJ, 
([2]). Properties of thermodynamics potential were also 
attributed to this (configurational) free energy, and it 
possesses the property of mininrality as to a variation of 
concentration of structural defects. It is easy to see that 
this property issues from the fact that the free energy, in 
a logarithmic scale and with an opposite sign, coincides 
within a constant with PDF ([7|) 



F,v = -kBT\n{Zf), 



(20) 



and as a consequence the minima of the free energy au- 
tomatically coincide with the maxima of PDF (compare 
curves 1 and 2 of Fig. 4), which determine the most prob- 
able, that is, equilibrium states. In other words, the free 
energy Fcv is a quantity reciprocal to PDF, expressed in 
energy units on a logarithmic scale. 

Energies F and Fcv to some extent similar have im- 
portant principal distinctions. Namely, it is known that 
the internal and free energies, related by the Legendre 
transformations, depend on the arguments of its own 
U = U{T, V), F = F{S, V), we name them symbolically 
the eigen-arguments of these potentials. This property 
and the respective differential relationships, allow to in- 
terpret these energies as thermodynamics potentials. If 
we compare the energy of vacancy uy to the tempera- 
ture T, then in this context, it is necessary to confront 
the number of vacancies ny with entropy S. It would 
seem the configurational free energy Fcv must be a func- 
tion of vacancy energy uy, however, from relationships 
([T]), ^ and (ITO)) it is follows that it is a function of the 
number or concentration of vacancies ny ■ It can be con- 
cluded that the free energy Fcv is not a thermodynamics 
potential in strict sense. 

To construct some thermodynamics potential of the 
free energy, note that in the case of vacancies according 
to (1) and (2), the configurational entropy Scv and the 
number of vacancies ny are in a mutual one-valued de- 
pendence with each other, that is, the number of vacan- 
cies nv can always be used as an independent argument 
instead of the configurational entropy Scv- In this case, 
executing transformation of Legendre type with the pair 
of the conjugated thermodynamic variables uy and ny, 
we get the alternative type of the free energy (curve 4 of 
Fig. HD 



F,^ 



It is easy to show 



nv 



U - Uynv- 



dF,^ 



duy 



(21) 



(22) 



Relationships (IT5|) and (1^^ are the joint pair of equa- 
tions for the internal energy U and the modified config- 
uration free energy Fcv on one side, and for the vacancy 
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FIG. 4. Basic set of thermodynamic functions for quadratic 
approximation: 1 - probability distribution function /; 2 - 
classic configurational free energy Fc; 3 - internal energy U; 
4 - a tangent to U at the equilibrium point; 5 - modified 
configurational free energy Fc; 6 - effective internal energy 
U; 7 - effective configurational free energy Fc- The right- 
hand scale is for PDF. Straight line corresponds to equilibrium 
state. 



concentration ny and vacancy energy uy on the other 
side. It can be seen that the vacancy concentration is an 
eigen-argument of the internal energy, and the energy of 
defects is an eigen-argument of the modified configura- 
tional free energy. In the quadratic approximation this 
dependence can be written in an explicit form (see curve 
4 of fig. H 



Fcv — U() 



1 



2uvi 



{uvQ -uvY 



(23) 



As the energy necessary for the formation of a new va- 
cancy is lower, in the presence of other defects, than in 
a vacancy-free crystal, the quadratic correction in P^ 
has negative sign. Note that expression (|12p is suitable 
for both the equilibrium and non-equilibrium state. In 
this approximation the internal energy is a convex func- 
tion of the vacancy number, having a maximum at point 
ny = nvmax, as is shown in the insert of Fig. |4l In that 
approximation the modified configurational free energy 
is a concave function with a minimum at point ny — 
(or uv = Mo)- 

With relationships P^ and (|22l) it is easy to show that 
the stationary state corresponds to neither the maximum 
of the internal energy U nor to the minimum of the free 
energy Fcv, as the stationary state is at point ny — nye, 
where uve and ny,, satisfy equations 



UVe 



dU 
dn 



Ve 



^0, 



dK 



nve = 



cV 



duve 



+ 0. (24) 



It is known from thermodynamics that the product 
TScv in (17), entered in determination of the canonical 
free energy Fcv, makes sense of the bound energy lost 
for the production of work. On the other side, the total 




n(%) 



FIG. 5. "Bound energy" of a solid with vacancies: curve 1 
PDF; curve 2 bound energy in form TScV', curve 3 bound 
energy in form uvny- 



energy of defects uynv, in the main part, physically is 
the energy which is also lost for work production. Only 
a small part of it remains for work production. Then we 
can conclude that these two energies must be equal, at 
least, in the region of the equilibrium state (Fig. [5]). 



TScv ~ uynv- 



(25) 



The energies are equivalent as we have deal, actually, 
with energy of the same origin, but written in different 
forms. The quantities, included in this relationship, are 
exactly enough defined, and energy wyny, is to be un- 
derstood as the total energy stocked in the defects. Tem- 
perature T is the temperature of thermomotion, and it 
is directly unconnected with a defect subsystem. Ap- 
proximate equality (23) reflects a dynamic thermal equi- 
librium between the thermal subsystem, as a dynamic 
motion of atoms, and the static subsystem of defects. 
The configurational entropy Scv is an interface between 
these subsystems. It is also interesting to note that both 
temperature T and energy of vacancy uy are specific en- 
ergetic characteristics of each subsystem. In addition, 
both the configurational entropy Scv, and defect concen- 
tration Tiy possess property of additiveness, that along 
with their mutual one- valued accordance ([T]), © can be 
foundation for their interchangeability in general thermo- 
dynamic relationships. 

Summarizing the above said, we can apply similar rela- 
tionships and reasoning to interstitial atoms. The ques- 
tion of their application to other types of defects (disloca- 
tions, grain boundaries etc) remains open, as there are no 
simple and reliable statistical expressions for them, anal- 
ogous to those presented in part II. At the same time, it 
is obvious that from the thermodynamics point of view 
all structural defects are characterized in identical way by 
surplus energy, spent for their formation, and they can 
be described by similar evolution equations, but written, 
for the present, at phenomenological level. 



D. Kinetic (evolution) equations 

If the system has deviated from the equihbrium (sta- 
tionary) state it will tend to return back to that state at 
a speed the higher the larger the deviation. Evolutional 
equations for a solid with point defects, in this case, are 

drim _ , , dU _ 

a, — -"-Trim I a ^me): 

at orim 



dUm 

dt 



Tlu^i 



dUr. 



'^mej ' 



(26) 



Hereafter m = V,i 

The form of kinetic equations (|26p is symmetric with 
respect to the use of internal and configurational free en- 
ergy. If the equilibrium state is closer to the maximum 
of internal energy, then, in view of solution stability, the 
upper sign is chosen (convex function). If it is near the 
minimum of internal energy, the lower sign is taken (con- 
cave function). Both variants of the kinetic equations 
are equivalent and their application is a matter of conve- 
nience. 

At the same time, for a number of reasons, more prefer- 
able are evolution equations expressed in terms of the in- 
ternal energy. First, concentration (or number) of defects 
is a more easily measurable variable, second, the inter- 
nal energy is a base quantity in thermodynamics (the first 
law of thermodynamics) , it is universal both for the equi- 
librium and non-equilibrium state and, finally, internal 
energy is the analogue of Hamiltonian. Setting the inter- 
nal energy, as a function of state variables we specify the 
system and fully determine it from the thermodynamics 
point of view. 

If we assume that the equilibrium energy of defect 
uv — uve (or Ui = Uie) and the number of defects 
ny = nve (or rii — riie) are slowly changing during an 
external load, then their product can be introduced un- 
der differentiation sign in (j26p . It permits us to specify a 
new kind (shifted) of internal and free (effective) energies 
(see curves 5 and 6 in Fig. |4]): 



U U UjyiQTljyi^ r Q^ 
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Then equations (P^ are simplified a little 
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(28) 



The original potentials U and Fcm are connected by 
means of transformation of the Legendre type (PT|) . The 
shifted potentials U and Fcm are connected by means of 
transformation 



^ cm — ^ '^m'^m i 



^Tl -\- U'fYiTl^i 



(29) 



which differs from the Legendre-type transformation by 
anticommutation bracket [wm"-m] = Ume^m +i*m'^me- 

For the shifted potentials the stationary point coin- 
cides with maximum of U (curve 5 of Fig. 2]) and with 




FIG. 6. Internal (1) and effective internal energy (2, 3). 
nve(i) and nve(2) are vacancy concentration for the first and 
second stable equilibrium state, Uf, is an inflection point of the 
internal energy (1). For calculations, parameters were chosen 
close to parameters of copper uvo = l{eV), uyi = 0.022(61^), 
UV2 = 0.00017(6^). 



minimum of F^m (curve 6 in Fig. |4]) for the upper signs 
in Eq. ((28|) . Thus, U is an effective thermodynamic po- 
tential, for which the tendency of the original part of 
the internal energy to a minimum is completely compen- 
sated by the entropic factor. Free energy Fcm tends to 
a minimum, but it tends to a minimum in the space of 
eigen-argument Um, while the free energy Fc tends to a 
minimum in the space of "non-eigen-argument" Um- 

In the presence of only one equilibrium state n,„e, Ume, 
the effective internal energy is determined simply. In 
the presence of two such states transformation ([77]) can 
done for each one. In this case, we get two graphs of the 
effective internal energy (Fig. |6l). 

As the first state is in the region of convexity of the 
internal energy (curve 1), the graph of its respective ef- 
fective internal energy has a maximum (curve 2), and in 
equations (PS|1 , (^5]) the upper signs are chosen. The sec- 
ond state is in region of concavity of the internal energy, 
and the graph corresponding to its effective internal en- 
ergy has a minimum (curve 3), and in equations (j26p . 
(PS)) lower signs are chosen. Obviously, the evolution of 
the system from non-equilibrium states to the nearest 
equilibrium states till the inflection point on the graph 
of internal energy should be described by evolution equa- 
tions (j28p . taken with upper signs, and past at that point, 
vice versa, with lower ones. 



III. DISCUSSION AND PROSPECTS 

Above it has been shown, that with the increase in the 
temperature of a solid some new state with high concen- 
tration of vacancies is formed, which in accordance with 
the hypothesis of Ya. I. Frenkel can be interpreted as the 
liquid state. At the same time, according to the Born hy- 



pothesis a distinctive property of liquid is the absence of 
shear rigidity or equahty to the zero of the shear modulus 
of material. Shear rigidity displays activation character, 
and it is in inverse dependence with mobility of point 
defects. 

To study the issue let us consider a solid having only 
vacancies, by setting the internal energy in the form 

u{nv) =uo + uvqUv - -uviuy + -Mv2"y, (30) 

where Uq , uyi {I — 1,2,3) are some coefficients depen- 
dent on temperature T and elastic deformation e^ , as on 
control parameters. For the sake of convenience and fur- 
ther standardization we have passed from the number of 
vacancies to their concentrations ny {ny — ny/N), and 
in a similar way to specific internal energy u ~ U/N. 
Then 

^A(4)^ + ^(e^,)^ 



liO 



2 
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UVQ 



^vo + ^^i^iV+fii^tjy 



(31) 



where A and /i are Lame moduli for a vacancy-free solid, 
A, p, are corrections to the elastic moduli due to exis- 
tence of vacancies (defect of a modulus) , e? is the first 
invariant of the tensor of deformations, {^fj)'^ = ^tj^'ji- 
It is assumed that contribution of invariants is the most 
substantial at the lower degrees of presentation (1501) . 

Expression for the effective shear modulus within the 
framework of this model is 



fief = fJ,~ fiynv- 



(32) 

The sign "minus" is chosen provided, that the effec- 
tive modulus in material with vacancies is less than in 
the vacancy-free one (see Fig. 2.3 in [49]). The first 
term fj, is shear modulus in a vacancy-free crystal. It 
can be interpreted as a local (microscopic) modulus of 
the lattice in regions far from vacancies. The effective 
(macroscopic) modulus ^e/ is a quantity averaged over 
volume containing a lot of vacancies. It is remarkable, 
that at a phase transit point at melting, namely, the ef- 
fective macroscopic modulus goes to zero exactly, while 
the microscopic modulus does not convert to zero |39]. 

A few types of the temperature dependence for the 
shear modulus are known. For example, there is a 
quadratic dependence ^{T) = /iO(l — T/Tn ) 13911. Other, 
more known semi-empirical relationship is |50l452l | 



m(T) = Aio - 



D 



exp 



To 
T 



(33) 
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where /ig is the shear modulus for T = OK, and D, Tq are 
material constants. The constant To can be treated as an 
activation energy Tq = Ea/ks, and the relationship p3p 
is presented as 

Ea 

Mo - (po + D) cxp {--r-^) 
m(T) = j^ ^^' . (34) 



^-^"P(^) 



According to Born principle, the temperature of melt- 
be defined on condition fi(Tm) — 
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A formal agreement between the theory presented in 
part II and relationship psp can be obtained by taking 
the temperature of phase transition for the temperature 
of melting T„j . At the same time, it is more logically to 
present the shear modulus as a two part one consisting of 
two terms, actually, shear modulus at zero temperature 
fiQ and an activation depressed term 



M(r) =Mo(l -cxp(- 



Eg 

knT 



))■ 



(36) 



However, in obedience to this expression, the shear 
modulus does not go to zero at limit temperature T^, 
but only ai T = oo. It is known in the theory of amor- 
phous materials, that activation energy is proportional to 
the shear modulus Ea ~ M 531 • It can be assumed that 



a similar supposition is just for crystalline solids too, at 
least, in the pre-melting state, for in this state its struc- 
ture becomes amorphous-like. Then relationship p6p can 
be rewritten as 



Me/ = Mo(l -exp(^ 



' knT 



))■ 



(37) 



It is taken into account that here we have the effective 
shear modulus. 

It is remarkable that this equation has a trivial zero 
solution /Lie/ = 0. This solution is important as it can 
serve for description of non-activation motion of liquid 
phase to which all solids pass as a result of melting. 

The graphs of temperature dependences of the shear 
modulus are shown in fig. [T] Curve 1 proper to semi- 
empirical theory p4l) . slowly changes at low tempera- 
tures and sharply diminishes to zero, passing to the re- 
gion of negative unphysical values. Curve 2, built for the 
effective shear modulus by formula ([37]), changes already 
at low temperatures and is more declivous. It smoothly 
tends to zero, not passing to the region of negative values. 
In addition, zero values of the effective shear modulus 
identically satisfy equation p7p , describing a liquid phase 
legalistically in the whole of temperature interval, how- 
ever, this solution is realized only after the phase tran- 
sition. Curve 3, which corresponds to the shear micro- 
modulus, does not go to zero at a "melting" point. 

The graphs of Fig. [7] correspond a continuous change 
with the temperature of the shear modulus in the solid 
phase. However, at the real melting-point because of the 
first-order phase transition the change can be sharper 
(see the arrow). Similar transitions are observed in the 
real materials (see, for example. Fig. 1 in Ref. |54|). 

At low temperatures, a basic contribution to shear de- 
formation is from dislocations in view of their greater spa- 
tial size and higher mobility. In obedience to the Cottrell 




T(K) 



FIG. 7. Temperature dependence of the shear modulus: curve 
1 semi-empirical shear modulus with (|34|) : curve 2 effective 
modulus with (|37|l : curve 3 shear micro- modulus. Arrow 
indicates the jump of the shear modulus to zero value at the 
real melting point. 



defect mobility must already be taken into account, as 
well as to take into account the influence of other types 
of defects on the process. 

Earlier in a phenomenological variant, the approach 
was applied for the description of processes of defect 
formation and work hardeniiigin metals treated by se- 
vere plastic deformation 30-33', auto- vibration transi- 
tions between the amorphous and nanocrystalline states 
in amorphous alloys [56] , processes of slipping and stick- 
slip in the super-thin lubricants [53, [5^. In particular, 
within the indicated models, to describe the stick-slip ef- 
fects the processes of transitions between solid and liquid 
states, that is, processes of melting and solidification, are 
examined. In prospect, the method can be applied for the 
description of kinetic processes of defect formation and 
forming the properties in different physical systems, in- 
cluding soft- matter [53, ISO] , complex protein compounds 
fH, etc. 



theory, mobility of dislocations is the higher, the wider 
their dislocation cores [53 • Mobility of vacancies at low 
temperatures is negligibly small. 

At more high temperatures the situation, however, 
changes cardinally. Vacancies due to the thermomo- 
tion become delocalized in space, formiiig the so-called 
extended defects of higher mobility |46l - |48| . The total 
energy barrier of a vacancy splits on a great number 
of barriers of lower height, which are considerably easy 
overcome due to thermal fluctuations. The number of 
variants of overcoming these barriers increases thus, here 
the vacancies have certain advantage over dislocations. 
The extended vacancy is a three-dimensional defect and 
the number of variants for "lateral" passing by potential 
barriers grows quadratically, while dislocation is a two- 
dimensional object and the number of variants in plane 
the perpendicular to line of dislocation increases only lin- 
early. Thus, in transition to the liquid state at point 
higher than the temperature of melting, a vacancy mech- 
anism becomes so perfect, that a solid loses resistance to 
the change of form fully. And the vacancy emptiness is 
uniformly distributed (dissolves) in the bulk of material 
so, that a single vacancy can not be practically identified 
and its presence is talked about only legalistically. 

For a more deep accordance it is necessary to take into 
account the fact of temperature dependence of the shear 
modulus in expansion of the internal energy (|5(I)) . (1511) . 
If we substitute relationships ^U^ ^^ to PDF dlTl), we 
get a more complicated dependence of PDF on tempera- 
ture in the presence of double exponents there. However, 
this does not at all affect the calculations of equilibrium 
values of point defects concentration in solids being in 
the free state, because in this state ef^ = 0. Due to it, it 
becomes possible to examine the phase transition and the 
change of elastic properties (shear modulus) during such 
transition separately. In the presence of external stresses 
and the resulting elastic deformations, the influence of 



IV. SHORT CONCLUSION 

In the article, within the framework of ideas of non- 
equilibrium evolutional thermodynamics, the vacancy 
theory of melting of solids is considered. With the ac- 
count of removal of degeneration of vacancy energies due 
to their interaction and approximation of this interaction 
by polynomial presentation of the internal energy it was 
possible to generalize Frenkel concepts and to describe 
the process of melting as jump-like Ist-order phase tran- 
sition. It is evident that the influence of interstitial atoms 
on the process of melting is not of principle, as the corre- 
sponding material curves lie higher than the same curves 
for vacancies and do not cross the fundamental curve in 
this temperature region. They can only a bit shift the 
proper stationary (equilibrium) solutions for vacancies. 
Influence of dislocations on the behavior of the system 
can play a role at relatively low temperatures. At tem- 
peratures close to the melting point, because of the size 
effect of increasing vacancy mobility, they become de- 
cisive in shear deformation of the material. The shear 
modulus here tends to zero. It is shown that apart from 
the tending to the zero of the shear modulus of the solid 
phase, the system has another state with the shear mod- 
ulus identically equal to zero, which is realized exception- 
ally for the liquid state. The transition to this state can 
be realized earlier, than the shear modulus reaches zero, 
that is, the real melting point will lie somewhat to the 
left from the point at which the shear modulus becomes 
zero. 

In addition, in the article on the example of solid 
with vacancies, a more general theory of non-equilibrium 
evolutional thermodynamics has been illustrated. It is 
shown that kinetic equations can be written not only in 
terms of the free energy, but also in terms of internal 
and modified free energy. The whole system of inter- 
connected non-equilibrium thermodynamic potentials is 
obtained, including, effective internal and modified free 
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energy. Intercommunication between classical and modi- 
fied free energy through the component of bound energy 
in the both forms is analyzed (fig. [5]). 
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evolution of shear modulus at melting. 
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